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ABSTRACT: 

Consideration was given to the problem of optimal control of parabolic equations. The existence solution of the 
considering optijnal control parabolic problem is proved. The gradient of the cost functional by the adjoint problem 
approach is obtained. Lipschitz continuity of the gradient is derived. 
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I. INTRODUCTION AND STATEMENT OF THE PROBLEM 

The optimal control problems governed by partial differential equations have developed very fast in the 
last 30 years, and it has brought a promising and vital researching domain to the subject of mathematics. The 
optimal control problems governed by partial differential equations concern many applications in physics, 
chemistry, biology, etc., such as materials design, crystal growth, temperature control, petroleum exploitation, 
and so on. The relative details can be seen in [1-4], and so on. The partial differential equations involved in 
these problems include elliptic equations, parabolic equations and hyperbolic equations [5-7]. The reference [1] 
considered the problem of optimization of the optimal control parabolic problem with control in boundary and 
right hand side of the equation. In the present paper, the control action in the initial and boundary conditions of 
the considering optimal control parabolic problem. 

Let the considered process be described in O. ^ = {{x ,t): x e {Q ,1) , 0 < t <T jby the following problem: 

u ^ = (x) u ^) ^ + V ^{x , t) , (x,f)Gf2,, , 

u{x,Q) = v,(x) , JC E {0,1), (1) 

Ji X'^,t) = 0 , - /I (/) u^{l,t)= Ic [u{l,t)- v^{t)] t e (0,r] 

where u {x,t) is the solution of the system (1), the constant ^: > 0 is called the convection coefficient or heat 
transfer coefficient , A{x) > 0 , A. {x) e L^[0,l] and V = V , x V ^ x V ^ is the set of admissible controls 
where V = {v = (v Jx, f ), (x), (0) : v Jx, 0 e J,v,(x)e (O, /), (0 e L^{0,T)} 

is a closed and convex subset in V c L, (Q ^ ) x L, [O, /] x [O, r ] . 

The problem of optimal control lies in determining admissible controls v e V minimizing together with the 
corresponding generalized solution of problem (1) the functional 

fA^) = P \ \u{x,T)- zix)fdx + a f^[v,(0- wit)Ydt (2) 

Jo Jo " 
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In (2), z{x) and w (f) are given functions, respectively, from [O,/] and I, [0,r], J is a fixed 
time and /? , a are given positive numbers. 

II. EXISTENCE SOLUTION OF THE OPTIMAL CONTROL PROBLEM 

In this section, we give the definition of the weak solution of the problem (1) and the existence solution 
of the optimal control problem (l)-(2). 

The weak solution of the problem (1) will be defined as the function u e L, (fJ , ) , which satisfies the 
following integral identity; 



j u {x,T ) {x,T)dx - j Vj (x) ^ (x,0)dx ~ [u{x,t) ^{x,t) ~ A{x) ^{x,t) ii ^{x,t)] dxdt 

- kf^[u{l,t)- v^{t)]c{l,t)dt = v^{x,t) C{x,t) dxdt ,V^(x,OgLJqJ (3) 

Evidently, under the above conditions with respect to the given data, the weak solution u e (f2 j- ) of 
the direct problem (1) exists and unique [8,9]. 

We define a solution of the optimal control problem (l)-(2), according to [10], as a solution of the 
minimization problem for the cost functional (v ) , given by (2): 

/.(^.)=inf /.(^) (4) 

veV 

Evidently, if (v, ) = 0, then the solution \\ E V is also a strict solution of the optimal control 
problem (l)-(2), since v, 6 V satisfies the functional equation _^ = z{x) , x e (O,/). Further, in 

the view of the weak solution theory for parabolic problems, one can prove that if the sequence {v'"'}c V 
weakly converges to the function v e V , then the sequence of traces [u[x,T ■,v^"'')] of corresponding solutions 
of problem (1) converges in (Q ^ ) - nonn to the solution {u(x, T; v)}, which means (v'"') (v) , 

as n— >oo. This means the functional (v ) is weakly continuous on V, hence due to the Weierstrass existence 

theorem [11] the set of solutions = {v e V : (v, ) = (/^ )^ = inf (v)} of the minimization problem 
(2) is not an empty set. 

ni. FRECHET DIFFERENTIABILITY OF THE COST FUNCTIONAL AND ITS GRADIENT 

The principal result in this section is Theorem 3.1. Its proof will be prepared by two lemmas. 
Let us consider the first variation 



' - - ' - ,2 



/„^')= + Av)- /Jv) =2/1 1^ [u{x,T-v)- z{x)] Au{x,T-v) dx +/? [Au(x,T;v)] 



dx 



+ a\ [Av (f)] dt +2a[ [v At) - w{t)\ At) dt (5) 
Jo Jo 

of the cost functional (2), where 

V + Av = {v„(x,r)+ Av^{x,t\v^{x)+ Av,(a:), v,(f)+ Av^(f)} e V 
Au{x , t; v) = u[x , t , V + Av) — u [x , t , v) <= L^[n. ^) . 
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Evidently the function Am = Au{x,t;v) is the solution of the following parabolic problem 



All ^ = {a{x) Au ^) ^ + AVg{x,t) 
Au{x,0) = A V| (x) 



, (x,f)s 
X e (O,/) 



(6) 



Am ^(0,f) = 0 
Lemma 3.1. 



A{i)Au ^{l,t)= k[Au{l,t)- Av^{t)] , ?G(0,r] 



Let V = {Vg{x,t),v^{x),v^{t)} , V + Av = {vg{x,t)+ Av^ {x,t\v^{x)+ Av^{x),v^{t)+ Av^{t)} e V be given 
elements. If u = u{x,t;v)e L^{n^^.) is the corresponding solution of the direct problem (1) and 
If (x,t;v) e L , (n , ) is the solution of the adjoint parabolic problem 



if{x,T) = ip [u{x,T;v)- z{x)] 



, {x,t) e Q. ^ 
, X s (O,/) 



j^_^(0,f)= 0 



A{l)if ^{l,t)= k y/{l,t) , t G (0,r] 



(7) 



then for all v £ V the following integral identity holds: 

2P\ [u{x,T \v) - z{x)] Au{x,T;v)dx = f i// {x ,0; v) Av {x) dx 
Jo Jo 

T 

+ I// (x,t-v) Av^^ {x,t) dx dt + k ^ I// {l,t;v) Av^{t) dt , VveV (8) 
Proof: Let us use the final condition at f = Tin (7) to transform the left-hand side of (8) as follows: 

2/? f [ u{x,T;v)- z{x)] Au(x,T;v) dx 
Jo 

= f if/{x,T;v) Au{x,T ■,v) dx 
Jo 

= jj [ I// {x,t;v) Au {x,t;v)]^dx dt 

= II [if ^[x , t; v) A u [x , t ; v) + If [x , t ; v) Au ^[x , t; v)]dx dt 

= II [- If ^{x , t;v)) Au{x , t;v) + if [x , t;v) Aw^(x,f;v))] dx dt 

nj 

+ II V ^ (x, /) If {x, f ; v) dx dt 

= { [- A {x ) If [x , t ; V ) A u [x , t ; V ) + If { X , t ; V ) A [x ) A u , { X , t ; V ) ] dt 
Jo ' 

+ II {x , t )if {x , t; v ) dx dt . 

Taking into account the boundary conditions in (7) and (8) for the functions y/(x,t; v) and A m (x , f ; v ) ; we obtain 
(8). 
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We will define the parabolic problem (7) as an adjoint problem, corresponding to the inverse problem (l)-(2). 
The parabolic equation (7) is a backward one, and due to the "final condition" at t = T it is a well-posed initial 
boundary-value problem under a reversal of time. 

Now we use the integral identity (8) on the right-hand side of formula (5) for the first variation of the cost 
functional (v) . Then we have 



A/^(v) = A Vg{x,t)i// {x,t;v)dx dt +k ^ i// (l , t; v) Av ^{t) dt H- j i// {x ,0; v) Av ^{x) 



dx 



+ 2a \ [vAt)-w{t)]AvAt)dt +a \ [Av At)] dt +p\ [Au{x ,T w)]' dx (9) 
Jo Jo Jo 

Taking into account the above definition of the scalar product in V and the definition of the Frechet- 
differential we need to transform the right-hand side of (9) into the following form: 



A/„(v) =(/;(v),Av)^ +«£ [Av,(0]"A [Au{x,T;v)\ dx (10) 

This formula provides further insight into the gradient of the functional (v) via the solution of the adjoint 
parabolic problem (7). Due to the definition of the Frechet-differential, we need to show that the last two terms 
on the right-hand side of (10) are of order O ( ||v ||^ ), with p > \ . 

The following result precisely shows an estimate for the last two terms in (1) in order O ( ||v ||^ ). 
Lemma 3.2. 

Let Au = All {x,t;v) G (fi j. ) be the solution of the parabolic problem (6) corresponding to a given 
V G y. Then the following estimate holds: 



f [Au{x,T;v)] dx < — Av , V Av e V (11) 

Jo g II llv 



where ||^^||y = I || l^^o l-*^' ^)| dx dt + ^ |aVj(x)| dx + ^ \Av^{t) \ dt I 
L", ° " J 

is the norm (f2 ^ ) - norm of the function A v e V , and the constants c^, , e >0 are defined as follows: 

, , , , f /' 2A: 1 
c = max , 1, = min X[x)> 0 , s = rem \ , \ (12) 

Proof. 

Multiplying both sides of the parabolic equation (10) by Au , integrating on Q , 

Am Am = — [a m ^ ]^ , (x ) A m ) Au ={x{x)Aii Au) - X{x) {Au Y and using the initial and boundary 

2 ' ' 

conditions we obtain: 

0= II Am [Am, - [X{x) Au ^) ^ - Av ^[x,t)]dx dt 
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= — j| [am ^ ]j dt dx + {a{x)Au ^ Au)^dx dt +jj A{x){Au dx dt +jj Vg{x,t)Au{x,t) dx dt 



2. 



1 ' 2 ^ 2 ^ 

— f [Au{x,T;v }] dx + k i [Au{l,t;v )] dt - k \ Au{l,t;v)Av{t)dt 
-7 Jo Jo Jo 



1 ' 



+ Jj"/l(x)(AM (x,/;v)) dx dt — v ^{x , t) A u (^x , t;v) dx dt — — J [AVj(x)] dx 



This implies the following energy identity 
1 ' 



J [a M (x , r ; V )] dx + /t J [a M (/ , f ; V )] dt + || 1 (x ) (A m J dx dt 



1 ' 



II Vj, (x, f) Am (x, f; v) d/ +fc| Am (/, f; v) A (f) ^ ~I )] '^■^ (13) 



for the solution A« = AM(x,f;v) of the parabolic problem (13). We use the e-inequality 

s a ~ 

a p < + ,\/ a , p e R , y s > 0 ,on the right-hand side integrals of this identity. Then we have 

2 2e 

T I I ^ 

jjvg{x,l)Au{x,l;v)dxdt +^| Au{l , t; v) Av ^{t) dt + — | [Ai'|(x)]dx 

s 0 1 1 ks ^ ■> 
< — II [Am (x, f ; v)]'ifedf + II [a v„ (x, ?)]'dxiif + 1 [Au{l,t\v)\ dt 

Q.J 2^Q^ ^ 

+ f [AvjIO] dt +— \ [Av,(x)] dx , Ve>0 (14) 

Further, we estimate the term [Am(x, ;)]^ by applying the Cauchy inequality, 



2 r ' 1 ' 

[AM(x,r)] = Au ^ , t; v)d — Au{l , t; v) ^ 

^ 2 

< 2[ I AM_j(^,f;v)d^ I + 2(AM(/,r;v))' 



2 2 

< 2/| [Am ^ (x, f; v)] dx + 2(Am (/, f; v)) 
Now integrate the both sides of this inequality on f2 j, 

T 

II [AM(x,f; v)]'dxA < 2/'|| [AM^(x,f; v)]'dxdf + 2/| [ Au{l , t; v)f dt (15) 
and use this estimate on the right-hand side of (14): 
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T 1 ' ^ 

^^Vg{x,t)Au{x,t;v)dxdt +kj Au{l , t; v) Av ^{t) dt + — | [Ai'j(.x:)]dx 



<£-/^ff[AM {x,t;v)]' dxdt +(£■/ + )f [Au(1 , t; v)]~ dt + [[ [Av {x , t)f dxdt 

JJ ' 2 2£ " 



k T , 1 ( 2 

+ f [Av^{t)Ydt + — f [Av,(x)] dx 

Is 2 •'° 

This inequality with (14) implies 

(l, - f /')[[ [Am .(j«:,f;v)]'dxdf + (yt - f Z - — ) f [Ak (/, f; v)]' A +— f \Au{x,T\v'^ 



1 , yt r , 1 ' ^ 

< ff [Av„(x,/)]"t/xc?f + f [Av^{t)]'dt +—\[Av(x)]dx (16) 

2e 2s ^0 ' 2 

ks 2 

Requiring the positivity of the terms k - si and A - si we get bound (12) for the parameter e >0. 

2 

With this parameter e >0, from estimate (16) finally we obtain: 

1' 2I 2 k 'T 2I' ^ 

— f [Au{x,T ■,v)] dx < ff [Av (a:,/)] dxdt + f [Av (/)] dt + —\ [Av ix)] dx 

2 2s 2s 2 

The required estimate (13) follows from this inequality by choosing the constant c^^ > 0 as in (16), which 
completes the proof. 

The lemmas 3.1, 3.2 imply that the last integral in (10) is bounded by the termO (||a v||^ ). Thus by the 
definition of Frechet-differential at y £ V 



= {fai^)'^''}^ +«£[Av,(0] dt +fiY[Au{x,T;v)] dx , 

We obtain the following theorem: - 
Theorem 3.1. 

Let conditions in the considered problem hold. Then the cost functional is Frechet-differentiable, (v ) £ 
C ' (V ) ■ Moreover, Frechet derivative at v e V of the cost functional (v) can be defined by the solution 
I// e W^" {Ci J ) of the adjoint problem (7) as follows: 

/J(v)= { V^{x,t-v), y^{x,();v)- k ^^(/,?;v)} (17) 
IV. LIPSCHITZ CONTINUITY OF THE GRADIENT 

Any gradient method for the minimization problem (4) requires an estimation of the iteration parameter a j. 
> 0 in the iteration process 
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= " -"JaV' I ' n=0,l,2. 



(18) 



where v e V is a given initial iteration. Choice of the parameter a j, defines various gradient methods [12], 
although in many situation estimations of this parameter is a difficult problem. However, in the case of Lipschitz 
continuity of the gradient / j (v ) the parameter a can be estimated via the Lipschitz constant as follows: 



0 < <y„ < a„ < 2/(L+ 2<yJ 

where 50, 51 > 0 are arbitrary parameters. 
Lemma 4.1. 

Let conditions of Theorem 3.1 hold. Then the functional (v ) is of Holder class c (v ) and 

||/;(v + Av)- f;{v)\l < l\\av\1. 



(19) 



(20) 



where 



||/;(v + Av)- /J(v)|[ = jj{Aiy{x,t;v)y-dxdt + i ' (A ^ (/, v))' 



dt 



+ ( — +— +— ) f {Aiy {x,0;v)) dx 
A, k 2 ^0 

and the Lipschitz constant L > 0 is defined via the parameters , e >0 in (16) as follows: 



L = 4 







/ k ' 


+ 


— + — 






k 2) 






x,t;v 4 


Av)- 


y/ {x,t; 



> 0 



(21) 



(22) 



Proof. 

The function A y/ {x,t;v) = 
backward parabolic problem: 

^A(i^, = (- /l(x)A(i^ ,(x,/)eQ, 
Aif/{x,T)= 2 Au{x,T;v) , xe(0,/) 

Aiy^XO,t)=0 , - A{l)Ay/ ^{l,t)= k Ay/{l,t) , fe(0,r] 



(23) 



Multiplying both sides of Eq. (23) by A y/ [x j\v) , integrating on O. ^ and using the initial and boundary 
conditions, as in the proof of Lemma 3.2, we can obtain the following energy identity: 

2 T 2 1 ' 2 

X{x)[Ay/ ^{xj;v)\ dxdt +k^ [a y/ {l ,t\v)\ dt + — j [A y/ {x ,Q\v)\ dx 



= 2 P \ [Au{x,T w)]^ dx 
Jo 

This identity implies the following two inequalities: 



(24) 



1 ', 



X,^^ [Ay/ ^{x,t\v)] dx dt + [Ay/ {x,Q;v)] dx < 2 /? j [Au{x,T w)]' dx 

T 7 1 ' 2 ' 1 

\ [Ay/ {l,t;v)\ dt + — \ [A y/ {x ,Q ; v)\ dx < 2 /? f [Au{x,T w)]' 
Jo 2 •'o 



dx 
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2 1 2 1 

Multiplying the first and the second inequality by and , correspondingly, summing up them, 

A, k 



and then using the inequality (15) we obtain: 



A. k 



f [Au{x,T w)]' 
Jo 



dx 



2 r I I 2 

II [a ij/ {x,t;v)] dxdt + ( + —) \ [a ly {x ,0 ; v)] dx <4fi 

JJ ^ ^ Jo 

Let us estimate now the second integral on the right-hand side of (21) by the same term. From the energy 
identity (24) we can also conclude 



f [Ay/ (l,t;v)ydt +— f [Aiy (x,0;v)] dx < 2 P k \ [Au{x,T w)]' 
Jo 2 •'o 



dx 



This, with the above estimate, implies 

2 ^ T 2 I k I 2 

ff [a [// (x,f;v)] dxdt + k' \ [A {y (/,f;v)] dt + ( + — + — ) f [a ^ {x,Q ; v)] dx 

•'° A, k 2 •'<' 



< 4/? 



^/^ Ik 



A k 2 



■ + — \ \ [Au{x,T \ v)\ dx 
2 ) 

Using this in (21) and taking into account Lemma 3.2 we obtain (20) with the Lipschitz constant (22). 

V. CONCLUSION 

Consideration was given to the problem of optimal control of parabolic equations. The existence 
solution of the considering optimal control parabolic problem is proved. The gradient of the cost functional by 
the adjoint problem approach is obtained. Lipschitz continuity of the gradient is derived. 
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